Light scattering by an oscillating dipole in a focused beam 
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The interaction between a focused beam and a single classical oscillating dipole or a two-level 
system located at the focal spot is investigated. In particular, the ratio of the scattered to incident 
power is studied in terms of the oscillator's scattering cross section and the effective focal area. Debye 
diffraction integrals are applied to calculate it and results are reported for a directional dipolar wave. 
Multipole expansion of the incident beam is then considered and the equivalence between this and 
the Debye diffraction approach is discussed. Finally, the phase change of the electric field upon the 
Q\ ■ interaction with a single oscillator is studied. 



PACS numbers: 42.50.Ct,03.65.Nk,32.50.+d,32.i 
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INTRODUCTION 

The realization of quantum networks and repeaters for quantum information science crucially depends on an efficient 
interface between photons and single quantum systems Strong interaction has been achieved by coupling single 
emitters to optical resonators 0, U 0] anc ^ nas a ^ so been predicted for emitters located in waveguides where the light 
^p" 1 . is tightly confined in the transverse dimensions El ^_ J la • In free space, the question is to what extent photons 
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12| . Recent experiments have demonstrated that light focused on 



may interact with a single oscillating dipole 



single ions, molecules, or quantum dots may be attenuated in transmission by a few percents 0,HQ0Q In 
q-( a recent theoretical study we have shown that a focused light beam can be perfectly reflected by a single oscillating 



dipole located at the focal spot [18| . In this paper we investigate in more detail the cases where the beam is a focused 
t-H ■ plane wave and a directional dipole wave. We discuss the equivalence between the multipole expansion and the Debye 
diffraction approaches. Moreover, we compute the phase shift induced on the electric field by the dipole and find that 
| a few degrees are easily obtainable using realistic focusing parameters and off-resonance excitation. 

The strength of the interaction between a beam and an oscillator can be expressed by /C, the ratio of the scattered 
to incident power. K. can also be given as the ratio of two independent quantities, the scattering cross section a and 
» the inverse of an effective focal area 

O. K= P~ = A' () 

. £_| \ a and A depend exclusively on the oscillator and focusing setup properties, respectively. For a classical oscillator and 
' a two-level system (TLS) the cross section reads 
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where we assumed that there is no damping other than by radiation, ctq = 3A /(2tt) depends solely on the wavelength 



of the transition 20( ■ r results from radiation reaction for the classical oscillator 20( , while A represents the Einstein 
coefficient of spontaneous decay 21] . A denotes the detuning from resonance and Q is the Rabi frequency of the TLS, 
which imposes saturation effects at stronger incident light. 

For a point-like oscillator the scattered power depends solely on the field strength at the position of the oscillator. 
Accounting for the electric nature of the interaction, the effective focal area can be given as the ratio of the power 



transmitted through the focal plane (FP) and the electric energy density at the focal spot 18 1 
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FIG. 1: a) The incident light propagating from left to right is focused onto an oscillating dipole located at the focal spot in 
vacuo. GRS: Gaussian reference sphere, a: entrance-aperture radius, a: entrance half angle, /3: collection half angle, /: focal 
length, O: focal spot, b) Transmittance as a function of the laser detuning is displayed for a focused plane wave (FPW) with 
a = f3 = 7r/3 and for a directional dipole wave (p x ) with a = j3 = n/2. 



where S z denotes the z component of the Poynting vector in the FP and W c i(0) is the electric energy density at 
the focal spot O 22| . The integration is taken over the FP. The second equality holds for circular symmetry of the 
incident field strength with respect of the z axis; a condition which, for instance, is obeyed by a focused plane wave 
(FPW) but not by a directional dipole wave. Figure la) describes an ideal lens that projects the incident field onto 
the Gaussian reference sphere (GRS), which represents the locus of equal phase of the incoming converging and also 
for the outgoing diverging mode. Because the lens is assumed to be in the far-field region, the fields are tangential on 
the GRS. 



DEBYE DIFFRACTION 



An established approach of calculating the field in the focal area is provided by the Debye diffraction integrals. 
This approach was initiated by Debye using Green's theorem 23] and was extended by Wolf using the method of 
stationary phase 24j. For an incident plane wave the method was extensively applied by Richards and Wolf [2z 



These considerations led to the Debye diffraction integral for the electric and magnetic fields in the focal area 

Ae lfcr - S d£ , H(r)=--^-/ sxAe' 4r5 (ffi 
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where A denotes the vectorial angular spectrum of the incident wave, k is related to the wavelength by k = 2tt/X. r 
is the position in the focal area and s is the unit vector in the direction of the plane wave. The integration is carried 
out over the incident solid angle S; nc bounded by the semiaperture angle a. Calculations are presented in Fig. 2 for 
the case of a directional dipole wave indicated by p^. Such a wave is constructed by considering the emission pattern 
at the left hemisphere of the GRS of an electric dipole located at O and oriented along the x axis and by reversing 



the propagation direction 26[ . The angular spectra of the FPW 22], |27| and the p^ wave read 



fE V cos 9 (cos</> eg — sin^e^ 
A = { fEo (cos 9 cos (j> e$ — sin <j> e^) 
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In Figs. 2a) and 2b), the Poynting vector component S z and the electric energy density proportional to \E X \ 2 are 
displayed along the x and y axes in the FP. Note that the field components E y and E z are zero on these axes. Figure 
2c) shows areas of positive and negative values of S z in the FP, i.e. areas of forward and backward prop agation. The 



changes of direction are a signature of field vortices in the focal area, as reported for a FPW [25|, |28| . In Fig. 2d) 
the phase of the electric field relative to that of a plane wave is plotted for positions along the z axis. One notes a 
characteristic phase anomaly in the neighborhood of the focal spot associated with a phase jump of ir, which is also 
termed Gouy phase [29l. |30| . Furthermore there are oscillations, which do not vanish for increasing z displacements. 
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FIG. 2: a) The z component of the Poynting vector S z (full curve) and the electric energy density given as |-Ex| 2 (dashed 
curve) of a p x wave along the x axis in the FP and normalized to their respective values at x — 0. b) Same as a) but along 
the y axis, c) Contour plot of S z in the FP. Bright and black areas refer to the Poynting vector in the positive and negative 
z direction, respectively, d) Phase of the focused electric field E x on the z axis relative to that of a plane wave. Note the 
oscillatory behavior that does not vanish for large |z|-vahies. a — n/2 in all cases. 



This behavior is singular for propagation along the z axis, while for directions increasingly tilted away from the z axis 
the oscillations progressively die out at larger distances [23]. 

Using Eq. © we calculated A for four different cases, namely for the FPW, the p x wave, the dipolar wave with the 
generating dipole oriented along z axis, and for combined electric and magnetic generating dipoles directed along the 
x and y axis, respectively |l8|. Here we present, pars pro toto, the results for the electric field E(O) at the origin, the 
effective focal area A, and the scattering ratio K. for the ~p x wave 

E(O) - -i^°T(a) e x , - = ^F(a) , ^ = *F(a) , (6) 

3 A 3n 

where the subscript of ICq indicates that a = cr is assumed. The p x wave has the property that the three quantities 
depend in the same way on the semiaperture angle a through 

J- (a) — — (4 — 3 cos a — cos 3 a) . (7) 

As pointed out in Ref. 18], /Co reaches for a = ir/2 the maximum possible value of 2, which also establishes the 
maximum possible scattering ratio for a directional focused beam in free space. ICq > 1 indicates that the scattered 
power is larger than the incident power. However, this does not violate the energy conservation law because of 
destructive interference in the forward direction. Taking the interference into account the transmittance T, the ratio 
of the transmitted and incident power, is given by 

T = l-K = l-\^, (8) 

where TZ is the reflectance, the ratio of the back scattered to incident power. The factor of 1/2 in the second equality 
accounts for the fact that equal amount of scattering takes place in the forward and backward directions. Based on 
the procedure outlined in Ref. [3] we also determined the transmittance as a function of the semiaperture angle a 
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FIG. 3: Transmittance To of a p x wave as a function of the angles a and (3 as denned in Fig. la. The dashed curve indicates 
the edge along the geometrical shadow boundary a = ft. 



and semicollection angle (3 

%(a, /?) = 1 — — (4 — 3 cos a — cos 3 a) (4 + 3 cos (max{a, /?}) + cos 3 (max{a, /?})) , (9) 
16 

where the subscript to 7o indicates that a = cto is assumed. Examples of T as a function of the detuning are presented 
in Fig. lb) for the FPW and for the p^ wave. Figure 3 displays a rapid decrease of Tg with increasing a and an edge 
along the geometrical shadow boundary a — f3, as for the FPW [3] . As shown in Eq. (J9j> , % is invariant with respect 
to j3 for (3 < a, while for (3 > a, T increases with p. Contrary to the FPW, 7o(a, (3) decreases monotonously with 
increasing a and reaches the value of zero at a — ir/2. 



MULTIPOLE EXPANSION 



Another approach convenient for the description of focused fields is given by a multipole expansion [27Ll31l.l32l.l33l.l3 



Adopting the notation of Bohren and Huffman [32j we write for the electric field in the most general form 

i 

E(r) = ^2 ( B e,m,<M e , mi <(r) + A eim ^N e>m ^(r) + e -> o) , 

l m=0 



(10) 



m £ (r) and Ne m ^(r) are real valued and denote complete sets of magnetic and electric multipoles. Be m t 



where M 

and Ae m £ are the corresponding coefficients. Assuming a linearly polarized field in front of the incident lens and 
aligning the x axes to the incident field polarization the expansion in Eq. IjlOp can be restricted to M 0j i^(r) and 
N ej i^(r) multipoles for the electric field and to M^i^r) and N 0) i j f(r) for the magnetic field, respectively. 

The calculation of the coefficients requires some attention j^ . The source- free field mode may be considered as a 
sum of the converging incoming and diverging outgoing field. Because the outgoing mode is purely a consequence of 
the incoming mode, only the latter is needed for a unique determination of the coefficients. This concept was applied 
for instance by Sheppard and Torok [27| , where the multipoles of the expansion were associated with spherical Hankel 
functions. The direct expansion in terms of multipoles for the source-free field is also possible. However, in this case 
the converging field at the entrance and the diverging field at the exit of the GRS have to be taken into account. 
For this purpose the field symmetry on the GRS has to be considered 36|, |37I | , which can be derived from the Debye 



scattering integrals in Eqs. Q when assuming positions diametral with respect to the origin 



Ae~ ikr s dS = 
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FIG. 4: Multipole expansion coefficients for the non-focused plane wave (PW) and the p x wave. \A e ,i,t\/Eo = \B ,i,e\/Eo 
(dots) are presented for the PW and \A £t i : e\/(Eofk) (dots) and \B B> \^\/ (Eofk) (circles) for the p x wave, respectively. For the 
latter a semiaperture angle of a — tt/2 was assumed. 



which means that the field is antihcrmitian for diametral positions on the GRS. The corresponding relationship of 
the field's phase (p reads 



y>{x, y, z) = —tp{—x, —y, —z) — tt, mod 2n 



(12) 



The phase shift of — tt demonstrates the phase anomaly in the neighborhood of the focal spot (see Fig. 2d) and it is 
equal to the Gouy phase acquired when the beam traverses the focus. 

Here we follow the approach of Borghi 38| and Borghi et al. 39( and expand the angular spectrum A of the 
incident field in surface vector harmonics and substitute the expansion into Eqs. (4). This procedure assures that the 
Debye-diffraction and multipole-expansion method are literally the same. We write for the expansion 

(-*) 



A = 



2k 



(13) 



where M Q - 



and fi emi i{9,<f>) are real valued and complete sets of vectorial surface harmonics, which are 



independent on the radial variable r. They are related to the spherical vector harmonics Y and Z by 

. / 2^+1 ^ 1/2 



Y 1 



\2Tii(e 

f 21+1 x 1/2 
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where Zj = s x Y^. Because of the completeness and orthogonality of the basis functions, the coefficients are given 
by 

2£+ 1 
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where the prefactors result from normalization of the basis functions and from accounting of the Whittaker type 
of transformation, which provides a relationship between surface vector harmonics and multipoles 4ll. |42|. For the 
magnetic multipoles this relationship reads 

(-»)' 



4tt 
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and analogously for the electric multipoles 
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FIG. 5: Angular spectra A of the incident field according to Eq. ((5| and approximations according to Eq. (|13|) are given by 
dashed and full lines, respectively, for a — tt/2. (a) FPW. (b) Dependence along the x and y direction, indicated by <j> = and 
<f> = it/2, respectively, for the p x wave. The inset indicates the angle 6 as it is considered for this figure. The summation of 
terms is truncated to f max = 8 and 16 in (a) and (b), respectively. 



We finally write the fields in terms of multipoles 

oo 

E(r) = ^(B oa ,*M , u (r)+A e , M N e , M (r)) , 

. oo 

H(r) = —J2( B o,U N o,M+A e ,i,tM et i,i(r)) , 
1=1 

and for completeness we also present expressions for the surface vector harmonics and multipoles 
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where ji = ji(kr) is the spherical Bessel function. Si is a related to jt, and iri and Tg follow from the Legendre 
polynomials Pj as 



S,= 



d(krj e (kr)) 
d(fcr) 



P/(cos( 
sin. 9 



dP/(cos( 

dT~ 



(19) 



In Fig. 4 we depict the coefficients and j5 ,i,£ of the p x wave for a = n/2, and compare them with the 

coefficients of a non-focused plane wave (PW) 32|. We note that A e ^j differs from zero for even I except for I = 1 
while the -B ,i,f coefficients differ from zero exclusively for odd I. The fact that coefficients with I > 1 do not vanish 
for the Pa; wave is somewhat surprising. However, these are required to maintain the propagation characteristics of a 
directional wave and to guarantee power conservation throughout the space on the basis of a source-free focused field. 
Figure 5 displays the quality of the expansion for the FPW and p x wave when the number of terms is truncated. It 
is apparent that quite a few terms are required for a decent reproduction of the angular spectra for a — n/2 and even 
more terms are required for a < it/ 2. 

Of special interest is the property that all multipoles are zero at the origin except the electric dipole mode N^ 1 ] x . 
Therefore, in cases where the field at the origin only is relevant, the analysis can be simplified by decomposing the 
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FIG. 6: Phase shift $ of Eq. (|26|) in units of degrees as a function of the laser detuning. Results are plotted for the FPW and 
for two different semiaperture angles a, as indicated. 



incident field into dipolar and nondipolar modes. This concept was introduced by van Enk 33j and was considered 



in Refs. 
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18l |. N e l l at the origin and in the far-field region reads [3 



N 
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e,l,l 



3 fa;, 



o , 



N (3) 
e,l,l 



i{kr--n/2) 

- — ^ (cos # cos </> e# — sin^e^) , kr ^> 1, z>0, 



(20) 



where the superscripts (1) and (3) denote the source- free field with the spherical Bessel function and the outgoing 
mode with the spherical Hankel function, respectively. The concept of field decomposition into dipolar and nondipolar 
components becomes particularly useful when the scattered field E sca is considered, which in the case of a classical 
dipole oriented along the x axis is given by 

a ikr 



E S ca(r) = -^inc(O)^^ ~— 



Inserting the expression 



(cos 9 cos 4> eg — sin <f> e^) , 



E inc (0) = 4,i,iNS, i(O) • e x , 



kr > 1 



(21) 



(22) 



into Eq. (|2"Tj) , it is fairly easy to see that at resonance the dipole component of the incident field is exactly canceled 
by the scattered field in the forward direction. Therefore, the outgoing field is given by 



E„ 



E;, 
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(23) 



2A + ir kr 

where, apart from Ej nc , the coefficient ^4 e ,i,i is the only variable that depends on the focusing specifications. For the 
FPW and the p^ wave, ^4 e ,i,i is given analytically as 



4i,i = 



j^fkE Q (8 - cos 3/2 a (5 + 3 cos a)) , FPW 
fkE (4 — 3 cos a — cos 3 a) , p x . 



TO 
.1 



(24) 



For a Pa; wave with a ~ 7r/2, A e ^ i/(fkE ) = 1 as shown in Fig. 4. Furthermore, by inserting Eqs. (j2TJ|) and (|24p into 
Eq. ([2^)1 one sees that the electric field at the origin is the same as in Eq. ©, confirming that the Debye-diffraction 
and the multipole-expansion methods yield identical results. 



PHASE SHIFT BY SCATTERING 

With the above considerations it is also easy to calculate the phase shift $ imposed on the beam by a single oscillator 
at the focal spot. The phase shift $ is defined by 

$ = arg(E out -E* nc ) , (25) 
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where E* nc is introduced as a reference field for the detection of the phase shift. Making use of Eqs. (|23| and (|24|) 
and assuming that a detector is positioned on the z axis, we find 

»T ttt (8 -cos 3 / 2 a (5 + 3 cos a)) \ .FPW 

<l> = ,u " I 1 - 2AT7T X V 3 ^ (26) 

2A + z7 l (4 _ 3 C os a -cos 3 a) J , P x , 

where an extra negative sign is introduced to account for the Gouy phase shift of tt imposed on the incident field. $ in 
Fig. 6 shows a typical dispersive of behavior. It amounts to 5-15 degrees at the extremal points for semiaperture angles 
accessible in experiments. The extremal points are located at approximately A/T = 1/2 and $ decays only slowly 
with increasing detuning. We expect that integration over a collection solid angle would not change substantially the 
picture gained from Eq. (f2"u) because of the coinciding phase fronts of the incident and scattered field. 



CONCLUSIONS 



We studied the scattering of a FPW and a p x wave by a single oscillator, with emphasis on the equivalence between 
the Debye diffraction and multipole expansion approaches. We systematically applied the concept of the GRS as 
the locus of equal phases in the forward and backward direction and paid special attention to the calculation of the 
multipole expansion coefficients on the basis of source-free fields. We thus derived an analytical expression 7o(a, /3) 
for the transmittance of a p^ wave. We finally demonstrated that a considerable phase shift of a few degrees is 
imposed on the light beam by a single oscillator at a detuning significantly larger than the linewidth. This property, 
for instance, might be exploited for the non-resonant detection of single emitters. 
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